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Abstract: We perform the Hamiltonian analysis of some form of the non- linear massive 
gravity action that is formulated in the Stiickelberg formalism. Following seminal analy- 
sis performed in arXiv:1203.5283 [hep-th] we find that this theory possesses one primary 
constraint which could eliminate one additional mode in this theory. We performed the 
explicit Hamiltonian analysis of two dimensional non-linear massive gravity and we found 
that this is theory free from the ghosts. 
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1. Introduction 

One of the most challenging problem is to find consistent formulation of massive gravity. 
The first attempt for construction of this theory is dated to the year 1939 when Fierz and 
Pauli formulated its version of linear massive gravity [Q] ^. However it is very non-trivial 
task to find a consistent non-linear generalization of given theory and it remains as an 
intriguing theoretical problem. It is also important to stress that recent discovery of dark 
energy and associated cosmological constant problem has prompted investigations in the 
long distance modifications of General Relativity , for review, see 

Returning to the theories of massive gravity we should mention that these theories 
suffer from the problem of the ghost instability, for very nice review, see |^. Since the 
General Relativity is completely constrained system there are four constraint equations 
along the four general coordinate transformations that enable to eliminate four of the six 
propagating modes of the metric, where the propagating mode corresponds to a pair of 
conjugate variables. As a result the number of physical degrees of freedom is equal to two 
which corresponds to the massless graviton degrees of freedom. On the other hand in case 
of the massive gravity the diffeomorphism invariance is lost and hence the theory contains 
six propagating degrees of freedom which only five correspond to the physical polarizations 
of the massive graviton while the additional mode is ghost. 

It is natural to ask the question whether it is possible to construct theory of massive 
gravity where one of the constraint equation and associated secondary constraint eliminates 
the propagating scalar mode. It is remarkable that linear Fierz-Pauli theory does not suffer 
from the presence of such a ghost. On the other hand it was shown by Boulware and Deser 
||5| that ghosts generically reappear at the non-linear level. On the other hand it was shown 

^For review, see M. 
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recently by de Rham and Gobadadze in [^j that it is possible to find such a formulation of 
the massive gravity which is ghost free in the decoupling limit. Then it was shown in Q 
that this action that was written in the perturbative form can be resumed into fully non- 
linear actions ^. The general analysis of the constraints of given theory has been performed 
in [Q]. It was argued there that it is possible to perform such a redefinition of the shift 
function so that the resulting theory still contains the Hamiltonian constraint. Then it 
was argued that the presence of this constraint allows to eliminate the scalar mode and 
hence the resulting theory is the ghost free massive gravity. However this analysis was 
questioned in |31] where it was argued that it is possible that this constraint is the second 
class constraint so that the phase space of given theory would be odd dimensional. On 
the other hand in the recent paper [^] very nice analysis of the Hamiltonian formulation 
of the most general gauge fixed non-linear massive gravity actions was performed with an 
important conclusion that the Hamiltonian constraints has zero Poisson brackets. Then 
the requirement of the preservation of this constraint during the time evolution of the 
system implies an additional constraint. As a result given theory has the right number 
of constraints for the construction of non-linear massive gravity without additional scalar 
mode ^. 

All these results suggest that the gauge fixed form of the non-linear massive grav- 
ity actions could be ghost free theory. On the other hand the manifest diffeomorphism 
invariance is lacking and one would like to confirm the same result in the gauge invari- 
ant formulation of the massive gravity action using the Stiickelberg fields. In fact, it was 
argued in |Q that for some special cases such a theory possesses an additional primary 
constraint whose presence implies such a constraint structure of given theory that could 
eliminate one additional scalar mode. However the limitation of this analysis is that it was 
not performed for the general metric so that one can still ask the question whether the 
elimination of the ghost mode occurs in general case as well. In fact, there is a well known 
example of the theory which seems to be consistent around some background while it is 
pathological for general background which is the Fierz-Paull theory Q which is ghost free 
around the flat background while it contains ghosts in general background [Q. This paper 
is the first step in the Hamiltonian analysis of Stiickelberg formulation of the non-linear 
massive gravity with general metric. More precisely, we consider one particular class of the 
non- linear massive action that was presented in ^. This action is sufficiently simple to 
be able to find the explicit relation between time derivatives of scalar fields and conjugate 
momenta while it possesses all interesting properties of the non-linear massive gravity as 
was shown in |32, ^ ^. 

We find the Hamiltonian form of given action and determine primary and secondary 
constraints of given theory. Then we show that due to the structure of the non-linear mas- 

^For related works, see 0, 0, 0, ^ |l| |l], ||, |l|, || g ||, ||, ^J|, |o[ 0. 

^Alternative arguments for the existence of an additional constraints were given in even if the 
Hamiltonian analysis was not complete and the minimal non-linear massive gravity action was considered 
only. 

*We should stress that the bimetric non-linear massive gravity is also manifestly covariant form of the 
non-linear massive gravity action | po| , In fact, it was argued in these papers that these manifestly 

diffeomorphism invariant massive gravities are ghost free as well. 
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sive gravity action this theory possesses one primary constraint that is a consequence of 
the fact that det = det{di(p'^dj(j)^ g'^) = 0) as was firstly shown in |2|] ^ . This result 



has a crucial consequence for the structure of the theory. On the one hand this constraint 
could provide a mechanism for the ehmination of an additional scalar mode however on the 
other hand the condition det V^^ = makes the calculation of the algebra of the Hamil- 
tonian constraints very difficult. In fact, we were not able to calculate this algebra in the 
full generality with exception of the two dimensional case where however two dimensional 
gravity is trivial. More precisely, the Hamiltonian treatment of two dimensional non-linear 
massive gravity shows that there are no physical degrees of freedom left and this result 



coincides with the analysis performed recently in pQ , 16 1. 

The structure of this paper is as follows. In the next section (|2|) we review the main 
properties of non-linear massive gravity and rewrite it into more tractable form that is 
suitable for the Hamiltonian analysis which will be performed in section @. In section @ 
we perform the Hamiltonian analysis of two dimensional non-linear massive gravity action. 
In conclusion (^) we outline our results and suggest possible extension of this work. Finally 
in the Appendix we perform the calculation of the Poisson brackets of the constraints 
when we do not impose the condition detF^^ = 0. Of course this is not the case of 
the non-linear massive gravity action but we include this appendix in order to show the 
complexity of given calculation even if we should again stress that this should be considered 
as a toy model calculation. 

2. Non-Linear Massive Gravity 

Our goal is to study non-linear massive gravity action in the version that appears in 

S = M^I d4xv^(4)ii(g) - ^M^m^ J d^x^/^gU{g-^ f) . (2.1) 

Note that by definition g^'^ and f^y transform under general diffeomorphism transforma- 
tions x'^ = x'^{x) as 

-S-Wf^l^. /;,(.') =/^w|^£,. (2.2) 

Now the requirement that the combination ^~^/ has to be diffeomorphism invariant im- 
plies that the potential U has to contain the trace over space-time indices. Further, it 
is convenient to parameterize the tensor /^jy using four scalar fields cj)^ and some fixed 
auxiliary metric f^uifp) so that 

= d^^^d.cl)'' fABicp) , (2.3) 



where the metric /ab is invariant under diffeomorphism transformation x'^ = x^{x') which 
however transforms as a tensor under reparametrizations of (f)^. In what follows we consider 
Iab = VAB, where tjab = diag(-l, 1, 1, 1). 



^ We should stress that in the previous versions of given paper the constraint det V ^ — was not taken 
into account and hence wrong conclusions were reached. 
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The fundamental ingredient of the non-linear massive gravity is the potential term. 
The most general forms of this potential were derived in 1^, ^. Let us consider the minimal 
form of the potential introduced in 

Uig,H) = -A({lCf- 




n>l I n>2 



(2.4) 



where we now have 



n=l ^ l\ I 



(2.5) 



and where {H'"')u = Ha^^H^^ . . . Hu"~^ ■ Using this explicit form of the potential we present 
an important observation related to the fact that the potential has to be defined using the 
trace over curved space-time indices which is a consequence of the requirement of the dif- 
feomorphism invariance of the non-linear massive gravity action written in the Stiickelberg 
field. Then we observe 



B 



= {5\ - A^s)i6^^ - A^^) = n^n^A , 
{HX = H^^Ki . . . i^^- = ni^nli . . . n^'r' ■ 

These observation implies that the potential can be written in the form 

\n>l ) n>2 

= -4(< >2 - < >) = -4{ic\f + 4ic\ic\ , 

where we defined 



(2.6) 



(2.7) 



B-^B-\J^B~-^B- y^-°) 

We mean that the trace over Lorentz indices is more convenient for the Hamiltonian treat- 
ment since we can easily implement the ADM decomposition of the space time metric. 
Explicitly, we use 3 + 1 notation [p4[ ^ and write the four dimensional metric components 



'For review, see ISq]. 
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as 



goo = -N'^ + Nig'^Nj , goi = Ni , gij = gij , 
^00 ^oi^^ ^ij ^ i, N'N^ 



(2.9) 



Note also that 4— dimensional scalar curvature has following decomposition 

= KijG'^'^'Kki + ^^^R , (2.10) 
where ^^^R is three-dimensional spatial curvature, Kij is extrinsic curvature defined as 

Kij = ^{dtgij - ViNj - VjNi) , (2.11) 

where Vj is covariant derivative built from the metric components gij. Note also that 
is de Witt metric defined as 

g^JM ^ }_^g^kgJl ^ gUgjk^ _ g^Jgkl g,,^^ ^ 1 ^^.^^^.^ ^ ^^^^ _ 1^^^.^^^ (2.12) 

Finally note that in ( ^.lOD we omitted terms proportional to the covariant derivatives which 
induce the boundary terms that vanish for suitable chosen boundary conditions. Using this 
notation we have 

A^B = -V„</.^V„</.B + g'^d^(t)^djct>B , V„0^ = ^{dtct^^ - N%<P^) . (2.13) 

We would like to stress that there is an another issue with the construction of the Hamilto- 
nian formalism. To see this note that the general potential term contains the matrix \/A 



that can be defined as power series y = X^„c„(A")^ with appropriate coefficients 
Cn- Then the variation of this expression is equal to 



SJA^J, = c^5A\ + ci{6A^c^'^B + aI-^A^b) + . . . (2.14) 



c 



that clearly cannot be written as 6yA^^ = ^JA"^ ^(\/A) . However the situation 

improves when we consider the potential term that depends on ^J~A^ since then we have 
7 



5^A\ = ci5Ai + 2c^6A%A\ + 2,c^6A\A''cA''A + ... = Ua^ ((VA)-!)""^ . (2.20) 



^Note that this resuh is consistent with the definition of the square root of the matrix given in ^] 

(Va)^s(Va)^c=A-^c (2.15) 

since 

S{^)% + (%/A)1,5(v^)^c((^)"')'d = <5(A)l3((%/A)-i)^o . (2.16) 

Taking the trace of given equation we immediately obtain ( 2.20 ). Note also that due to the matrix nature 
of objects A and B the following relation is not valid 

VaB^VaVB (2.17) 
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For that reason we restrict ourselves to the Hamiltonian analysis of the non-linear massive 
gravity action with the potential term 

U = 4TY\/A . (2.21) 
In fact, the Hamiltonian analysis of the gauge fixed form of this non- linear massive gravity 



was firstly performed in |31| and then in |32, |33[. It was shown in these two papers that 



this theory possesses interesting properties that allow to eliminate one additional scalar 
mode. Then our goal is to analyze given theory written in the manifest diffeomorphism 
invariant form and try to identify possible additional constraints that could eliminate one 
scalar mode. Explicitly, we would like to perform the Hamiltonian analysis of the following 
action 

S = Ml (f^dtN^ iKijO'^'^^Kki + (^)i? - m^Tr^j . (2.22) 
3. Hamiltonian Analysis 



In this section we perform the Hamiltonian analysis of the action ( 2.22| ). For the General 



Relativity part of the action the procedure is standard. Explicitly, the momenta conjugate 
to A^, Ni are the primary constraints of the theory 

7rAr(x) PS , 7r^(x) ^ (3.1) 

while the Hamiltonian takes the form 



(3.2) 



where vr'-? are momenta conjugate to gij with following non-zero Poisson brackets 

{5., (x), 7r'='(y)} = i (5f 4 + 6l5^) 5(x - y) . (3.3) 
Finally note that 7rAr,-7r* have following Poisson brackets with N,Ni 

{iV(x),^jv(y)} = 5(x - y) , {iV.(x),^-'(y)} = 5l5(^ - y) . (3.4) 



unless A and B commute. On the other hand since obviously A and A ^ commute the equation ( |2.17| ) 
gives 

VaVA-1 = I (2.18) 

which implies following important relation 

' = Va^ . (2.19) 
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Now we proceed to the Hamiltonian analysis of the scalar field part of the action. Note 
that in 3 + 1 formalism the matrix A"^ takes the form 

A^B = -V„0^V„,/.B + g'^dict>^d,ct>B = K^B + Vi , (3.5) 

where 

K^B = -^n^^^n<pB , Kab = VAcK% = KbA , 



Then the conjugate momenta pA are equal to 
PA- —V9^{^ )c- 



(3.6) 



(3.7) 



Note that using the symmetry of Aab = A ba we can write (^]^) in simpler form 

PA = M2m2^(A-i/2)^sV„0^ . (3.8) 
Using this expression we derive following relation 

-^^^PAPB = (A-V2)^^(v„c/>^V„</.^)(A-i/2)^^ = 



(A-V2)^^(yc;D_ACi>)(A-i/2)^^ 



(3.9) 



which implies 

HaB = (A-V2)^^yCD(A-l/2)^^ ^ (3_-LQ) 

where we introduced the matrix Hab defined as 

^AB = \irviPAPB + VAB ■ (3.11) 

gm^Mp 

Note that when we multiply ( 3.10| ) by V from the right we obtain (we use matrix notation) 

UV = (A"i/V)(A~i/V) (3.12) 

which implies 

A^^/V = VUV . (3.13) 

This relation will be important below. The important point of our analysis is that V^^ 
has the rank 3 as was firstly explicitly stressed in [p9| ]. In fact, if we introduce the 4x3 
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matrix = di(f>^ and its transpose matrix (W'^Y = di<f>^ which is 3 x 4 matrix we can 
write 

V^^ = Wjg'^{W^)j^ . (3.14) 

Then since W^,g^^ have the rank 3 we obtain that V^^ has the rank 3 as well. As a 
result det V = 0. In other words V is not invertible matrix. This fact has an fundamental 
consequence for the Hamiltonian structure of given theory. On the other hand we can 
multiply ( 3.10| ) by V from the left so that 

VU = {VA-^/^){VA-^/^) 

(3.15) 

that now implies 

VVU = VA-^/'^ . (3.16) 
With the help of these results it is easy to determine corresponding Hamiltonian 

= dt^^pA - Use = M>2^iVy^i?(A-l/2)^^ + N'PA^^C|>^ = 



(3.17) 



using ( 3.13| ) and using an obvious relation Tr\/VTT = Tr^/llV. With the help of these 
results we find the final form of the Hamiltonian formulation of the action ( 2.221 ) 



H = J d^^iNUr + N'Tii) , (3.18) 

where 

-Ht = n^^ + , -Hi = nf^ + nf , (3.19) 

where the explicit forms of these terms is given in (|3.2| ) and in ( ^.17 ). Then note that 
the requirement of the preservation of the primary constraints Tr^v w , vr* ~ implies an 
existence of the secondary constraints 

^t(x) w , ni{x) « . (3.20) 

The crucial question is the existence of additional constraints in the theory. The existence 
of these constrains would be important for the elimination of an additional scalar mode. 
We should expect that this constraint is the primary constraint between momenta and 
coordinates of the scalar fields and the question is whether we can find such a relation. For 
example, we can try to calculate det H. Using 

det (vAB + ^^.PAPb) = det ^ac det [s^^ + ^^,P^Pb) = 



det 7] exp Tr In ( H —j — tP^Pb] = det r] exp In ( 1 H —j — jPAP 



1 

gMjm'i 

(3.21) 



1 + .^4 4 PAP 
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and hence from (|3.1CI| ) we obtain 

1 A _ detV 

^ gM^m^^^^ ~ "ditA ' 

(3.22) 

We see that the upper equation imphes the primary constraint on condition when det V = 0. 
Then (|]2|) imphes the primary constraint 

C ■■ -T^PAP^ + 1 = 0. (3.23) 

The fact that det V^^ = is however crucial for the calculation of the algebra of the 
constraints. Unfortunately we are not able to perform this calculation in case of the 
four dimensional case due to the absence of the inverse matrix V~^. We demonstrate 
the complexity of such a calculation in the appendix when we abandon the condition 
det V^^ = 0. However we are able to perform this calculation in case of two dimensional 
massive gravity which we perform in the next section. 

4. Two Dimensional Massive Gravity 

Two dimensional massive gravity is exceptional also from the fact that the gravity is trivial 
so that the dynamical content is hidden in the scalar degrees of freedom only. In order 
to determine the physical number of degrees of freedom we have to find an appropriate 
structure of constraints. To proceed we have to determine the form of the Hamiltonian 
constraint ( |3.17 ). An explicit calculation gives 



(4.1) 



where we introduced following notation 



gii=uj , det 5 = cj , g^^ = - , di = d . (4.2) 



As a result the extended Hamiltonian takes the form 

H = j dxiNTir + N^^s + u'^ir'^ + u^ttn + uitt^ + TC) , 



(4.3) 

where we included the primary constraint vr"^ ~ into the definition of the spatial diffeo- 
morphism constraint. Note the this theory possesses following four primary constraints 

7r;v~0, vri^O, tt^ ^ , C= \ + 1 = . (4.4) 
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First of all the requirement of the preservation of the primary constraints 7rAr,7r^ implies 
following secondary ones 

-Hr ~ , ns^O . (4.5) 

We see that Tir ~ 0, Tis ~ are the secondary constraints. Then we have to check 
the consistency of all primary and secondary constraints during the time evolution of the 
system. The time evolution of the constraints tttviVt^ does not generate new conditions 
while the requirement of the preservation of the constraint tt^ implies 

-"<«» + ^'^^"-")!-rA.-rA = o, 



where we used following Poisson bracket 



(4.6) 



1 



{C(x),7r"(y)} = — pAP^(x),5(x - y) = A(x)5(x - y) (4.7) 

and we used the fact that we have to determine the evolution on the constraint surface. 
On the other hand the time evolution of the constraint C is equal to 

dtC = {C, H] ^ {C, Tt{N)] + j cLxvi^ {C, 7r'^(x)} = {C, Tr(A^)} + u^A = , (4.8) 

where we used the fact that 

{Ts{N^),C] = -N^dC (4.9) 
which vanishes on the constraint surface and where {C, TTi^)} is equal to 

{T.(iV). C} = 28 (EMpaiL) ?! + 28 (n'-P, ?± . 

V wv A / w V V A / 

(4.10) 

Note that this Poisson bracket vanishes on the constraint surface Us ~ , vr"^ ~ as well. 



As a result we find that (4.6) and ([4. 8)) determine that Lagrange multipliers u'^, T are equal 
to zero. 

Finally we have to determine the algebra of constraints Ht, Ti-s- As usual we calculate 
the Poisson brackets of their smeared forms and we find 



{TT(M),TT(iV)} = \j dxdy^^ {A{x),A{y)} ^^^^ 



A{x) TaM 



(4.11) 



To proceed we calculate following Poisson bracket 
{A(:r),A(y)} = 



-4.^{pAd^^){x) [d<P^{x)pA{y)dy6{x - y)- 
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- pA{x)dcl>^(,y)dJix - y)] -±-{pAd^^){y) + 

+ 4-i^(p^a0^)(a;) [{d<t>^{x)d<t>A{y)dy5{x - y)- 

- {d<l)%)dMx)dJ{x-y)\ -l-^{pAd<j>^){y) . 



(4.12) 



In the calculation of the Poisson bracket of the smeared Hamiltonian constraints we use 
the fact the Poisson bracket {A(x), A(y)} contains derivative of the delta functions that 
give non-zero contribution when they act on A'^ and M respectively. As a result we obtain 

{Tt{M),Tt{N)} = [ dx{MdN - NdM)^ X 
J Aa; 

= Ts (j^{NdM - MdN)^ + 2 J dx{NdM - MdN)uV'K'^ . 

(4.13) 

The algebra of constraints Tg takes the standard form 

{Ts{M'^),TsiN^)} = TsiM^dN'^ - N^dM^) . (4.14) 
Finally we determine the Poisson bracket 

{T5(M1),Tt(A^)} = j dx{Ts{M^),UT{x)]N = 

= - J dx{M^dHT + UrdX^) = TriM^dN) 

(4.15) 

using 

{Ts{M^),nTix)} = -M^&Ht - HtOM^ . (4.16) 

Finally we have to analyze the stability of the constraints T-Lt and Us- In fact, the 
constraint Ts{M^) is preserved due to the fact that the Poisson brackets between all 
constraints vanish on the constraint surface. On the other hand the time evolution of the 
constraint is given by the equation 

dtTriM) = {TTiM),H} « J dx(u'^ {TT(M),7r'"} + r{TT(M),C}) = 

(4.17) 

using the fact that u'^ and T vanish. 

At this stage we have finished the analysis of the time development of the constraints 
with following result. We have four first class constraints 

TTJV « , TT^ « 0, ~ , Til^O (4.18) 
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and the second class constraints 

tt'^ ^0 , C . (4.19) 

Using the second class constraints we can eliminate and w as functions of pA- Then the 
gauge fixing of the four first class constrains completely eliminate the physical degrees of 
freedom A^, vtat, A^^, tti and pA,4>^- 



5. Conclusion 

In this section we outline our results and suggest the possible extension of this work. We 
performed the Hamiltonian analysis of some particular model of non-linear massive gravity 



action written in the Stiickelberg picture [30|. We found corresponding Hamiltonian. Then 



following 1 29 1 where it was explicitly shown an existence of the primary and corresponding 
secondary constraint in the version of the non-linear massive gravity action presented in 
p3| ] we find corresponding primary constraint of the theory. Unfortunately due to the fact 
that this constraint is a consequence of the singularity of the matrix V^^ we were not able 
to determine algebra of all constraints and identify secondary constraints for the case of 
four dimensional non-linear massive gravity action. On the other hand we were able to 
complete the Hamiltonian analysis of two dimensional non-linear massive gravity where we 
showed that the algebra of the Hamiltonian and spatial diffeomorphism constraints is in 



agreement with the basic principles of geometrodynamics [36, 37, 38 1. We also identified 
an additional constraints and we show that there are no physical degrees of freedom with 
agreement with the analysis presented in [30, 



It is very unhappy that we were not able to finish the Hamiltonian analysis of four 
dimensional non-linear massive gravity theory due to the singular nature of the matrix 
V^^ especially in the light of the very nice proof of the existence of the primary and the 
secondary constraints that was performed in ||2^ in the case of the version of the non-linear 
massive gravity action presented in [^3|. To finish such an analysis is very desirable and it 
is the main goal of our future work. It would be also very nice to perform the Hamiltonian 
analysis of the general form of the non-linear massive gravity action with the Stiickelberg 
fields and we hope to return to this problem in future. 



A. Toy Model: Algebra of Constraints in Case of det V ^ 

In this appendix we perform the calculation of the algebra of constraints in case when we 
do not impose the condition det = even if these calculations are not directly related 
to the case of the non-linear massive gravity action studied in the main body of this paper. 
The goal of this appendix is to demonstrate the complexity of given analysis. To begin 
with note that the fact that det V^^ ^ implies the existence inverse matrix V^^ . In this 
case we find an important relation 

V~^^/Y^ = K-^l'^ = ^JmlV-^ (A.l) 

which will be useful when we calculate the algebra of constraints. 
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Let us consider the smeared form of the constraints ( |3.2C1| ) 



TriN) 



(A.2) 



The goal of this section is to determine Poisson brackets among these constraints. Note 
that in case of the General Relativity part of the constraints we have following Poisson 
brackets 



^gr{^)^A^ - y) - ^GR{y)^A^ - y) 



d 



n 



GR, 



d 



d 



(A.3) 



The calculation of the Poisson brackets that contains scalar phase space degrees of freedom 
is more involved. However it is easy to find the Poisson bracket between generators of 
spatial diffeomorphisms 



{?^f(x),7^f(y)} 



(X) ^^(X - y) - nr(y)l-S{x - y) 



(A.4) 



that together with the Poisson bracket on the third line in ( [A.3| ) implies following form of 
Poisson bracket between smeared form of the diffeomorphism constraints 

{Ts{N'),Ts{M^)} = TsiN^djM' - M^djN') . (A.5) 

In case of Hamiltonian constraint the situation is not so easy. Explicitly, from the definition 
of the Poisson bracket we find 

dVyNi^) W^(x),1^ff(y)}M(y) = 

A/ 



{T^"(7V),Tf (M)} = 



5pxiz) 



(A.6) 



To proceed further we use the relations 



6px{z) 



-^—{s^pg+paS^.)v^^ ( (ny)-V2^ ^ 



6{Vvu)^ciy 



lJ{VIif^{y) 
2 



B 



(x)5(x 



(z) 



((yn)-i/2)^^(y) 



1 



-W'dy.6{z - y)6'xdy,<P^iy)+g'^dy.cl>'^{y)dy,5piz - y)]UED ( (Vn) 



-1/2 



D 



G 



(y) • 

(A.7) 
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Then with the help of ( |A.7| ) we can determine the Poisson bracket (|A.6| ) . First of all these 
Poisson brackets contain the derivative of the delta functions. We perform integration by 
parts. Then the non-zero contribution arises from the derivative of the smeared functions 
N and M. As a result we obtain 

{Tf (A^),T^"(M)} = d^x{NdiM - MdiN)g'^ x 



AB 



= j d^^iNdiM - MdiN)g'^pA(yy/(UV)-Y''iy/(nV)-^U)Bcd, 

= j d^^iNdiM - MdiN)g'^pA{Vinvy^U)^Bdj(^^ = 
= T|^((iV5jM - MdjN)g^') , 



(A. 



where we used the symmetry of H and V and the relations 

\/(yn)-iy = v^J{Ilv)-^, nV(^n)-i = ^J{iiv)-m (a.9) 



that follow from ( A.l ) . If we combine this result with the Poisson brackets between smeared 
form of the General Relativity Hamiltonian constraints we find the final result 

{Tt{N), Tt(M)} = Ts{{NdjM - MdjN)g^') . (A.IO) 

Finally we calculate the Poisson bracket between T5(A^*) and T^(A^). Using 

{Tsm,pAi^)} = -N%pa{^) - diN'pAi^) , 

{T5(iV*),V9(x)} = -N'd,^ix)-m'^{x) , 
{Ts{N"),g^^{x)} = -N^dug'\x) + duN^g^^{x) + g'^dkN^{x) , 
{Ts{N'),Uab{^)} = -N'diUABi^) , 



(A.ll) 



With the help of these results it is easy to find 



(A.12) 



Collecting this result with the Poisson bracket between diffeomorphism constraint and 
Hamiltonian constraint of General Relativity we obtain 

{Ts{N'),Tt{N)} = TridkNN'^) . (A.13) 
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To conclude, we found that the theory with det V ^ possesses four constraints ?^t(x) ~ 
, Hi (x) ~ which are the first class constraints as follows from the Poisson brackets 
( |A.5| ),( |A.10| ) and (A. 13). We also showed that the crucial presumption for these calculations 



was the regularity of the matrix V^^ which of course is not the case of four dimensional 
non-linear massive gravity. 
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